Implementation of controlled phase shift gates and Collins version 
of Deutsch-Jozsa algorithm on a quadrupolar spin-7/2 nucleus 
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In this work Controlled phase shift gates are implemented on a qaudrupolar sys- 
tern, by using non-adiabatic geometric phases. A general procedure is given, for 



implementing controlled phase shift gates in an 'N' level system. The utility of such 
controlled phase shift gates, is demonstrated here by implementing 3-qubit Deutsch- 



> 

^ ' Jozsa algorithm on a 7/2 quadrupolar nucleus oriented in a liquid crystal matrix. 

o' 

§ ■ 1. INTRODUCTION 

o 

The use of quantum systems for information processing was first introduced by Benioff 
In 1985 Deutsch described quantum computers which exploit the superposition of multi 
particle states, thereby achieving massive parallelism 3|. Researchers have also studied 
the possibility of solvingcertain types of problems more efficiently than can be done on 
conventional computers j^, 0, Q . These theoretical possibilities have generated significant 
interest for experimental realization of quantum computers Several techniques are 

being exploited for quantum computing and quantum information processing, including 
nuclear magnetic resonance (NMR) 0, [9J. 

NMR has played a leading role for the practical demonstration of quantum gates and 
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algorithms. Most of the NMR Quantum information processing (QIP) experiments have 

nnnn 

utilized systems having indirect spin-spin couplings (scalar J couplings) [10j, [III, ]12|, 113] . 
Recently NMR QIP has been demonstrated in quadrupolar and dipolar coupled systems, 
obtained by orienting the molecules in liquid crystal media {14I In case of homo nuclear 
dipolar coupled systems, spins are often strongly coupled and hence can not be addressed 
individually lf|. However the 2 n eigen states are collectively treated as an n-qubit system 
Similarly for quad rup olar sys tems (s pin > 1/ 2), individual spins are treated as a 



multi-qubit system 
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281 ] . Resolved resonance 



lines provide access to the full Hilbert space of 2 n dimension |15l . 

A quadrupolar nucleus of spin I has (21+1) energy levels, which are equi-spaced in the 
high magnetic field and the 21 single quantum transitions are degenerate. In the presence 
of first order quadrupolar interaction this degeneracy is lifted and gives rise to 21 equi- 



spaced transitions 
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liquid crystals, such a situation is o 



be treated as an n-qubit system 



301 ] . For molecules partially oriented in anisotropic media, such as 
ten obtained [15[. If (21 + 1) = 2™, such a system can 
171 ] . The advantage of such systems and the dipolar 
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coupled systems, over the J-coupled systems is that the coupling values are one to two orders 
of magnitude larger, allowing shorter gate times or the use of shorter transition selective 
pulses. So far quadrupolar systems have been used for, quantum simulation, preparation 
of pseudo pure states, implementation of quantum gates and search algorithms |l8 
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. Recently Das et al have implemented Cleve version of 2- 
qubit DJ algorithm on a spin 7/2 nucleus j^. In all these cases the controlled gates are 
implemented by inverting populations between various levels, by using transition selective 
7r pulses. Recently it has been demonstrated that non-adiabatic geometric phases can be 
used for implementing quantum gates [33|, ]34j. Here we use non-adiabatic geometric phases 



to implement controlled phase shift gates and Collins version of DJ algorithm on a 3-qubit 
system obtained by eight eigen states of a spin-7/2 quadrupolar nucleus. 

The Hamiltonian of a quadrupolar nucleus partially oriented in liquid crystal matrix, in 
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the presence of large magnetic field B a , having a first order quadrupolar coupling, is given 

by, 

H = H Z + H Q = uj z + ^^9_S(3I 2 - I 2 ) = ujJ z + 2tt A (3/ 2 - J 2 ), (1) 

where uo Q = ~^B is the resonance frequency, S is the order parameter at the site of 
the nucleus, and e 2 qQ is the quadrupolar coupling. If the order parameter S is small, the 
effective quadrupolar coupling A can be of the order of a few kHz in spite of e 2 qQ being of 
the order of several MHz. Thus, it is possible to observe the satellite transitions due to first 
order quadrupolar coupling. 

A 50-50 mixture of Cesium-p ent adecafluoro-o ct anat e and D2O forms a lyotropic liquid 



crystal at room temperature 181 ] . The schematic energy level diagram of oriented spin-7/2 
nucleus is shown in Fig. flI]A.). The 8 energy levels are labeled as the basis states of a 
three qubit system. In the high field approximation, effective quadrupolar coupling (A) can 
be considered as a small perturbation to Zeeman field. Thus for population purposes the 
equilibrium density matrix can be considered to be proportional to H z (Fig. HA). Partially 
oriented 133 Cesium nucleus (1=7/2) gives rise to a well resolved 7 transitions spectrum at 
room temperatures ranging from 290 K to 315 K Fig (p]B). The effective quadrupolar coupling 
(A) changes with temperature, since the order parameter is a function of temperature. All 
the experiments have been carried out here at temperature 307 K, which gives A = 6856 
Hz. The equilibrium spectrum is obtained by applying a (7c/2) y pulse, is shown in Fig.(TjJ>). 
The integrated intensities are in the expected ratio 7:12:15:16:15:12:7, as determined by the 
transition matrix elements of \{Ix( y ))ij\ 2 (Appendix) [291 ] . 

2. CONTROLLED PHASE SHIFT GATES 



Non-adiabatic geometric phase in NMR was first verified by Suter et.al 32J. Non adiabatic 
geometric phases in NMR, were used to implement controlled phase shift gates, Deutsch- 
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Joz sa ( DJ) and Grover algorithms in weakly J-coupled and strongly dipolar coupled systems 
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A two level subspace (r,s) forms a fictitious spin- 1/2 subspace, hence the states |r) and \s) 
can be represented on a Bloch sphere (Fig. [2]) 34] . The two 7r pulses applied on a transition 



(r,s), cyclically rotates each of the states |r) and \s) 



34j |. The cyclic path followed by each 



34). The 



of the two states, makes an equal solid angle (Q) at the center of the Bloch sphere 
geometric phase acquired by the states |r) and \s) is given by e if ^ 2 and e~ %n l 2 respectively 
341 ]. where the phase difference of e™ between the states indicate that the two states 
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are traversed in the opposite direction |32|, |33J . It is possible to obtain various solid angles, 
hence various phases, by shifting the phases of the two ir pulses with respect to each other 



(Fig. [2]) [33j, |34j. It is shown that two 7r pulses with phases 9 and {6 + n + (p) respectively, 



applied on a transition (r,s), subtends a solid angle fl = 20, thus the phase acquired by the 
states |r) and \s) is e*^ and e~ 1 ^ respectively. 

However the internal Hamiltonian evolves during the application of a transition selective 
7r pulse, and gives rise to an additional phase known as dynamical phase. In order to observe 
only the geometric phase, one has to refocus the evolution of internal Hamiltonian. As shown 
in Eq. (pQ), the Hamiltonian consists of Larmor frequency term and quadrupolar coupling 
term. The Larmor frequency is identical to the frequency of central transition 'd'. Thus the 
receiver frequency is set at frequency of transition 'd' (Fig. [T]A.), hence there is no evolution 
due to the Larmor frequency term. The duration of the selective ir pulse is chosen in such 
a way that the phase acquired by the quadrupolar term makes a complete 2tc rotation, in 
other words e~ lHQtp = 1, where t p is duration of the selective ir pulse which is equal to 1.425 
ms. 

The phase information can be encoded in the coherences, hence in order to observe 
the phases the initial state should contain the coherences, which is obtained by apply- 
ing a non-selective (ir/2) y pulse on equilibrium state. Figure ([3^) shows the imple- 
mentation of controlled phase shift gate, represented by the diagonal matrix Ujl ^ = 
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diag[l, 1, —1, —1, 1, 1, 1, 1], obtained by applying two selective it pulses on transi- 
tion (3,4), with phases y and (y + -n + tt) = y. The two states |010) and 1 1 1) acquire tt 
phase shift so the transitions 'b' and 'd' are inverted, confirming the phase gate, since the 
phases of two tt pulses are same, it is possible to combine the two tt pulses in to a single 2tt 
pulse. Such pulses are used in section (3.1). 

In some quantum computing protocols (for example Grover search algorithm and quan- 
tum Fourier transform) one requires a phase gate such as (Table. 1), which creates 
the relative phase e 1 ^ at k th state, for example u£=[l, 1, 1, 1, 1, 1, 1, e % ^\. Such phase gates 
can be implemented by sandwiching various phase gates ut^ as shown in Table (1). 
indicates the geometric phase shift gate obtained by applying a pair of selective tt pulses on 
a transition (i,j) with phases y and (y + tt + </>), thus the diagonal elements corresponding to 
i and j states, acquire a phase shift e 1 ^ and e~ 1 ^ respectively. For example, u^-[ 8 2 ^=diag[e 1 ^^ 8 , 
e -^/ 8 ) i ) i ) i ) i ; i ; i] ; which can be obtained by applying two selective tt pulse on transition 
(1,2) with phases y and (y + tt + 0/8). 

The method of constructing controlled phase shift gates (Table 1), is explained here. 
For example, consider a system consisting of N energy levels and (N-l) single quantum 

transitions between the levels (1,2), (2,3), (N-1,N), as shown in Fig.(jl]). In this system, 

the controlled phase shift gate = diag[li, 1 2 , , (e l,? %, , ljv], can be implemented 

by sandwiching (N-l) phase shift gates, as, 



TT<f> _ ,.-*l N „- 2 */ JV _ _ o -(k-l)<t>/N (N-k).4>/N (N-(k+l)).4>/N _ 

u k — "(1,2) -"(2,3) • - u (k-l,k) ■ a (k,k+l) ■ u (k+l,k+2) ■ • U (N-1,N)\ Z J 

In Eq. d2J), by putting N=8 one can obtain any of the U% shown in Table (1). 

In Table (1), each of the operators, U%, requires seven pairs of selective it pulses, where 
each pair of pulses will act on a single transition. However the it pulses on unconnected 
transitions can be applied simultaneously. Hence one can simultaneously implement «(i,2), 
w (3,4)j w (5,6)j w (7,8) by using a pair of multi- frequency (MF) pulses, similarly W( 2 ,3), W(4,5), M( 6j7 ) 
can be implemented by using another pair of multi-frequency pulses. For example Iff can 
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be implemented by using seven geometric phase shift gates (Table 1), obtained by using two 
pairs of multi-frequency tt pulse, as described below. 



Jj4> _ .70/8 30/4 50/8 A/2 30/8 0/4 0/8 

U l ~ "(1,2) -"(2,3) -"(3,4) -"(4,5) -"(5,6) -"(6,7)- "(7,8) v / 

The unitary operators of phase shift gates correspond to diagonal matrices which com- 
mute with each other, hence Jjf of Eq. Q can be rewritten as, 



ut = {[«^.[«^.[«^.[4^]}.{[«%^[t.^].[<i]} 

— f\(>rr) a ' c ' e > 9 (ir) a ' c ' e ' 9 ]\ 
\[\ n )y "V'U (j/+7r+70/8),(j/+7r+50/8),(y+7r+30/8),(3/+7r+0/8) it' 

{[{TT)y '^•( 7r )( y +{- + 3^/4),(y+7r+0/2),(y+7r+0/4)]} (^) 

The 7r pulses acting on unconnected transitions, a, c, e and g (Eq. HJ), can be combined 
(added) in to a single pulse, since the operators of the subspaces corresponding to these 



transitions, 1%$ , ^(y) i ^iy) an< ^ ^x(y) 1^41' commu t e with each other. Similarly one can 
combine the (rr) pulses acting unconnected transitions b, d and f (Eq. @J. Figure ([3b) shows 
the spectrum obtained by a (ir/2) y pulse followed by the implementation of Uf(4> = tt/2) 
using four MF-7T pulses (Eq. HJ), the state 1 000) hence the transition 'a' acquires a 7r/2 phase 
shift, confirms U^ 2 . 



3. COLLINS VERSION OF DJ ALGORITHM 

Deutsch-Jozsa (DJ) algorithm determines in a single query, whether a given function is 



constant or balanced 



351 ] . A function is "constant" if it gives the same output for all 



inputs, and "balanced" if it gives one output for half the number of inputs and another for 
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Table (1). Unitary operators of controlled-^ phase shifted gates 





= 1,2, ,8) 


Uf = diagie*, 1, 1, 1, 1, 1, 1, 1] 


_,70/8 30/4 5^/8 0/2 30/8 0/4 0/8 


C/| = dm 5 [l,e^,l,l,l,l,l,l] 


_ ? -*/8 30/4 50/8 0/2 30/8 0/4 0/8 
- a (lj2) • «( 2> 3)- "(3,4) • "(4,5)- "(5,6) • "(6,7)- "(7,8) 


C/f = dmc/[l,l,e^, 1,1, 1,1,1] 


_ ,.-^/8 ? .-0/4 50/8 0/2 30/8 0/4 0/8 


L^ = dm 9 [l,l,l,e^,l,l,l,l] = 


_,,-^/8 ? -<A/4 -30/8 0/2 30/8 0/4 0/8 
- U (1>2) .«( 2 ,3) -"(3,4) -"(4,5) -"(5,6) -"(6,7) -"(7,8) 


^ = dm 9 [l,l,l,l,e^, 1,1,1] = 


-0/8 -0/4 -30/8 0/2 30/8 0/4 0/8 
- U (1>2) .M( 2 ,3) •'"(3,4) •'"(4,5)- M (5,6)-'"(6,7)-'"(7,8) 


^ = dm 5 [l,l,l,l,l,e-^,l,l] = 


- 7/-^ /8 7/"^ 4 ? /-^ /8 7/" 0/4 7/" 50/8 v/^ 4 «* /8 

- U (lj2) -U( 2) 3) -"(3,4) -"(4,5) '"(5,6) ' "(6,7) ' "(7,8) 


Uf = diag[l,l,l,l,l,l,e^,l] = 


-0/8 -0/4 -30/8 -0/2 -50/8 -60/8 0/8 
"(1,2) -"(2,3) '"(3,4) -"(4,5) '"(5,6) -"(6,7) ""(7,8) 


£/| = ^[1,1, 1,1, 1,1,1,6^] = 


-0/8 -0/4 -30/8 -0/2 -50/8 -30/4 -70/8 
"(1,2) '"(2,3) -"(3,4) '"(4,5) ""(5,6) '"(6,7) ""(7,8) 



the remaining half Classically for an n bit binary function, at least (2 n_1 + 1) queries are 
needed to determine whether the function is constant or balanced, whereas the DJ algorithm 
requires only a single query 6j. The Cleve version of the DJ algorithm requires an extra 
qubit (ancilla qubit) and uses controlled-not gates [35], whereas Collins version does not 



require the ancilla qubit but needs controlled phase shift gates 



36] 



For a three qubit system there are 2 constant and 70 (c^ 2 =c|) balanced functions [36l. 13 
The quantum circuit of the Collins version of the DJ algorithm is shown in Fig. ([5]). The 
algorithm starts with a pure state (pseudo pure state in NMR) 1 000) which is converted 

to a superposition state of eight basis states 1 000) , |001), ,|1H), by applying a pseudo 

Hadamard gate on all the three qubits. Thereafter a unitary operator U (oracle) is applied, 
followed by detection. Theoretically, after U one has to apply the Hadamard gate. In NMR, 
a Hadamard gate is replaced by a pseudo Hadamard gate, which is implemented by a (n/2)y 
pulse, and for detection another (ir/2)_ y pulse is needed. These two pulses cancel each other 
and the result of the algorithm is available immediately after U. 

The unitary operators (U) of the Oracle, are eight dimensional diagonal matrices. For 
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two constant functions, U's are given by, 



U c i = I (Unit matrix) and 

U c2 = diag[-l, -1, -1, -1, -1, -1, -1, -1] = (-l).J. (5) 

For balanced functions, we define the unitary operator as U(l,k,l, m), which means that 
each of the four diagonal elements l,k,l, and m are equal to -1, whereas the remaining four 
diagonal elements are 1. For example [7(1,2,3,4) = diag[— 1, — 1, — 1, — 1, 1, 1, 1, 1]. The 
unitary operators of 35 balanced functions can be written as, 



U(l,k,l,m) = [-|1)(1| - \k)(k\ - - \m)(m\].I, 
where k < I < m, k — 2, 3.., 6, I > k + 1, and m > I + 1 (6) 

The unitary operators of the remaining 35 balanced functions are identical to one of the 



operators of Eq. ([6]), up to an overall phase factor e m 371 ] . Here we implement the unitary 
operators corresponding to 11 of the 35 balanced functions (Eq. [6]) and one constant function 
(U cl ). 

3.1. Experimental Implementation 

(i) Preparation of pseudo pure state (PPS) 1 000) : A multi frequency pulse with six 
harmonics with appropriate amplitudes was applied. The six harmonics are the frequencies 
of the six leftmost transitions b,c,..., g of Fig. ([HB). Under the influence of this pulse, the 
populations of all the states except |000), start a collective population transfer among them 
as a linear chain of coupled oscillators 3|. Choosing the correct amplitudes of various 
harmonics and duration of the pulse, these six transitions can be simultaneously saturated 
with average population of these seven levels. The relative amplitudes of the six harmonics 
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b, c, and g are respectively given by 0.84, 0.93, 1, 1.03, 1.04 and 1.07 [27]] . The duration 
of the pulse is 2.05 ms. A gradient pulse was applied subsequently to kill (dephase) the 
coherences created during the process. The final populations were measured using a non- 
selective small-angle (5°) pulse. The small-angle pulse, converts the population differences 
between adjacent energy levels into single quantum transitions, within linear approximation. 
The spectrum of Fig ([6b) confirms the preparation of 1 000) pseudopure state. 

(iii) Coherent Superposition: After the creation of 1 000) PPS, a pseudo Hadamard 
gate (hard (n/2) y pulse) creates a superposition of 8 basis states. It may be noted that 
unlike in weakly coupled spins, the state created here is not in equal superposition of eigen- 
states, since the coefficients of various eigenstates are different. However, as is shown here, 
the created coherent superposition can be utilized for quantum parallelism, to distinguish 
different classes of functions. The state of the system after a ij/2) y pulse on 1 000) PPS, is 
given by, 



(e -i * J »).|000) 
1 



[|000> + V7|001> + >/2l|010) + V35|011> + «s/35|100) 



8V2 

+v / 21|101) + v / 7|110) + |lll)] 

= -^W|l) + V7\2) + V21\3) + V35\4:) + y/tt\5) + V21\6) + y/7\7) + |8>]. (7) 
8v 2 

The corresponding density matrix can be written as, 
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128 



|1> 
/ 1 



V7 



'21 
\/35 
y/35 

1 



|2) 



v/245 



V147 
7 

v/7 



13} 
21 



Vl47 



21 



V245 V735 



v/735 
21 



V147 
^21 



|4) 

\/35 
v/245 



v 7 ^ 



35 



35 



V245 
V35 



|5) 
35 



V245 
V735 



35 



35 



V735 V735 



V245 
V35 



|6) 
21 



VT47 
21 



21 



VU7 



1 7) 
7 

\/T47 



V245 V147 



V7S5 Vl47 



V147 



|8> 
1 \ 
V7 
TL 

/35 

/2l 



^7 



1) 

2) 

3) 

4) 

5)' 

6) 

7} 



where the elements within the boxes represent single quantum transitions a, b,.., g. Figure 
([6b ) shows the spectrum of where the intensities of transitions are obtained by modulus 
of the product of single quantum elements of Eq. [HJ, with corresponding matrix elements of 
4- 

(iv) Implementation of U: For constant function U = U c \ = I (Unit matrix), which 
requires no pulse, thus the final state if) c i is given by, 



\M = u, 



cl I 



(9) 



As mentioned in section (2), the controlled phase shift gate can be achieved by 

applying two (71") pulses with same phase on transition (i,j), which can be be combined in 
to a single (27r) pulse. For balance functions, U = U(l,k,l,m) (Eq. [6]), can be decomposed 
in to, 



U(l,k,l,m) = uj m .uj ltm) . 

— [ U Jl,2)---- U Jk~l,k)][ U Jl,l+l)- U Jl+l,l+2)--- U Jrn-l,m)] (^) 
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In Eq. ([TO]) . uj x ^ and uj t m % are implemented by applying two (2n) y pulses on transitions 
(l,k) and (l,m) respectively, if (l,k) and (l,m) does not correspond to single quantum tran- 
sitions, then they can be decomposed in to a series of single quantum transitions, as shown 
in the second equality of Eq. (flOl) . for example, 

[7(1,2,3,4) =u[ 1)2) .^ 3)4) = (27r) a .(27r) c , 

[7(1,4,5,8) = uJ 1A yUj 5 ^ = [M(l i 2)- M (2,3)-' U (3,4)]-[' U (5,6)-' U (6,7)-' U f7,8)] 

= [(27r) a .(27r) 6 .(27r) c ].[(27r) e .(27r) / .(27r) 9 ] 

= (27r) a ' c ' e ' 9 .(27r) b ' / , (11) 

where in the last equality of Eq. fill I) , the pulses acting on unconnected transitions (Fig. 
[Q are combined in to a single MF pulse. The required pulses for various U(l, k, I, m) of Eq. 
(J6]) are given in Table (2), the 2tt pulses acting on unconnected transitions, are combined in 
to a single multi-frequency (MF) 2tt pulse. The duration of each transition selective pulse, 
is set such that the evolution due to quadrupolar coupling, makes a complete 2tt rotation. 
It is to be noted that, in Eq. (TiTj) and Table (2), one can use any phase for (2tt) pulses, we 
have used 'y' phase for all the pulses. The final state \ipi } k,i,m) is given by 

\i(>i Mm ) = U(l,k,l,m)\il>). (12) 

It can be seen that l^i.fc^m) is same as \ip), except that the basis states |1), \k), \l) and 
|m) acquire a phase factor e ln (=-1). 

(V) Detection: The single quantum transitions (a, b, ..g) of the final states (Eq. |9][T2~]) 
are detected. Figure ([7]) shows the spectrum of \ip c i) and various \i[>i t k,i,m) ■ From the shape of 
the spectrum one can conclude, whether the final state represents a constant (or) a balanced 
function. For constant function {\ipci)) none of the peaks are inverted, whereas for balanced 
functions (\ipikim)) atleast one of the peaks is inverted. Furthermore the phases of transitions 
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Table (2). Unitary operators of balanced functions (Eq. and corresponding pulse sequences, 



consisting of 2tt pulses, where all the pulses are applied with phase 'y'. 



U(l,k,l,m) 


Pulse sequence 


U(l,k,l,m) 


Pulse sequence 


U(l,k,l,m) 


Pulse sequence 


U(l,2,3,4) 


(2vr) a > c 


U(l,2,3,5) 


(2vr) a ' c -(27r) rf 


U(l,2,3,6) 


(27r) a ' c ' e -(2vr) d 


U(l,2,3,7) 


(2^y,c,e_ ( 27r )^/ 


U(l,2,3,8) 


( 27r )a,c,e,g.( 27r )d,/ 


U(l,2,4,5) 


(27r) a ' d 


U(l,2,4,6) 


(2vr) a ' (i -(27r) e 


U(l,2,4,7) 


(27r)°' d -/-(27r) e 


17(1,2,4,8) 


(2vr) a ' d ^-(2vr) e ^ 


U(l,2,5,6) 


(27r) a ' e 


U(l,2,5,7) 


(2vr) a ' e -(2vr)/ 


U(l,2,5,8) 


(2 7r )«.e,s.(2 7r )/ 


U(l,2,6,7) 


(2vr) a ^ 


11(1,2,6,8) 


(2^) a '/-(27r)9 


U(l,2,7,8) 


(27r) a 'S 


U(l,3,4,5) 


(27r) a ' d -(2^) b 


U(l,3,4,6) 


(2vr) a ' d -(2^) b ' e 


U(l,3,4,7) 


(2^) a ' d '- f -(2 7 r) b ' e 


U(l,3,4,8) 


(2Tv) a ' d ' f -{2Tv) b ' e ' 9 


U(l,3,5,6) 


(27r) a ' e -(2vr) 6 


U(l,3,5,7) 


(2^) a ' e -(2vr) 6 '- f 


U(l,3,5,8) 


(27r)a>e,S-(27r) 6 >/ 


U(l,3,6,7) 


(27r) a '/-(2^) b 


U(l,3,6,8) 


(27r) a ^-(27r) b -» 


U(l,3,7,8) 


(2vr) a '3-(27r) 6 


U(l,4,5,6) 


(27r) a ' c ' e -(2vr) 6 


U(l,4,5,7) 


(27r) a - c ' e -(2vr) fe ^ 


U(l,4,5,8) 


(27r) a ' c ' e ' ff -(2 7 r) 6 > / 


U(l,4,6,7) 


(2 7r )°. c »/-(27r) 6 


U(l,4,6,8) 


( 27r )a,c,/_ ( 27 r) 6 '9 


U(l,4,7,8) 


(2^) a ' c 'f-(2vr) 6 


U(l,5,6,7) 


(2^) a ' c '- f -(2yr) M 


U(l,5,6,8) 


( 27r )a,c,/. ^M.S 


U(l,5,7,8) 


(27r) a ' c '9-(27r) M 


U(l,6,7,8) 


(27r) a ' c ' e ' ff -(27r) 6 > d 







confirm the final state iV'i^m)- F° r example in the spectrum of 1^1,2,3,4)5 transition d is 
negative since the phase difference between the states |4) and |5) is e l7T , while all others have 
zero phase difference, similarly one can confirm the other states. 

4. CONCLUSIONS 

Demonstration of quantum computing protocols on various NMR systems, is a promising 
research area for increasing number of qubits. In this work controlled phase shift gates are 
implemented on a oriented spin-7/2 nucleus, using non-adiabatic geometric phases obtained 
by using selective pulses on single quantum transitions. It is also demonstrated here that 
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in an N level system, one can implement a controlled phase shift gate by sandwiching 
various geometric phase shift gates, this method can also be applied to weakly as well as 
strongly coupled spin-1/2 systems. The number of selective pulses are reduced by using multi 
frequency (MF) pulses. Collins version of 3-qubit DJ algorithm is implemented, where the 
eight basis states are collectively treated as a three qubit system. The required controlled 
phase shift gates of the algorithm, are implemented by using MF-(27r) pulses. 



APPENDIX 



The angular momentum operators of spin-7/2 nucleus, are given by 29] 



I z \I,M) = M\I,M) 

1+ \I, M) = - M)(I + M + 1)| J, (M + 1)) 

/_ \I,M) = v/(J + M)(J-M + l)|/,(M-l)), (13) 

where 1=7/2 and M= -I, (-1+1),....,+! = -7/2, -5/2, -3/2, -1/2, 1/2, 3/2, 5/2, 7/2. 



1 1) 

2' 2/ 


11 5\ 
1 2' 2/ 


11 3\ 
1 2 ' 2 / 


11 1\ 

1 2 ' 2 / 


1 7 




11 n3\ 
1 2' 2 / 


11 =5\ 
1 2' 2 / 


1 2' 2 / 






A 

2 






















o\ 


1 7 
1 2 


Z ) 
21 





5 
2 






















1 7 
1 2 


^) 
21 








3 
2 



















1 7 
1 2 


1) 











1 

2 
















1 7 
1 2 


\) 














-1 

2 













1 7 
1 2' 


-1 



















-3 
2 








1 7 
1 2' 


-3 
2 






















-5 
2 





1 7 
1 2' 


-5 
2 


v° 






















1) 


1 7 
1 2' 


2 / 



(14) 



14 





/ „ 



ft 

V ' 
2 
















^ 




V I 
2 





f— 






















V3 





\/lE 
2 














I++I- _ 








VTE 

2 





2 











2 











2 





2 






















VTS 
2 





Vs 






















Vs 





2 




v° 

















2 





(15) 



( 


-is/7 
2 
















\ 


2 






























-iVIE 
2 




















2 





-2i 




















2i 





-iVTE 

2 




















iVTE 

2 





-iV3 




















iV3 





-i\/7 
2 


v° 

















is/l 
2 





(16) 



[1] P. Benioff, Quantum Mechanical Models of Turing Machines That Dissipate No Energy, Phys. 

Rev. Lett. 48, 1581 (1982). 
[2] D. Deutsch, Proc. R. Soc. London, Ser. A 400, 97 (1985). 

[3] D. Deutsch and R. Jozsa, Rapid solution of problems by quantum computation, Proc. R. Soc. 

Lond. A 439, 553 (1992). 
[4] P. W. Shor, Polynominal-time algorithms for prime factorization and discrete algorithms on 

quantum computer SI AM Rev. 41, 303-332 (1999). 



15 



[5] L.K. Grover, Quantum Mechanics helps in searching for a needle in haystack Phys. Rev. Lett. 
79, 325 (1997). 

[6] M.A. Nielsen and I.L. Chuang, "Quantum Computation and Quantum Information". Cam- 
bridge University Press, Cambridge, U.K. 2000. 

[7] M. Scholz, T. Aichele, S. Ramelow, and O. Benson, Deutsch-Jozsa Algorithm Using Triggered 
Single Photons from a Single Quantum Dot Phys. Rev. Lett. 96, 180501 (2006). 

[8] D.G. Cory, A. F. Fahmy, and T. F. Havel, Ensemble quantum computing by NMR spectroscopy, 
Proc. Natl. Acad. sci. USA 94, 1634 (1997). 

[9] N. Gershenfield and I. L. Chuang, Bulk spin-resonance quantum computation, Science 275, 
350 (1997). 

[10] I. L. chuang, L. M. K. Vanderspyen, X. Zhou, D. W. Leung, and S. Llyod, Nature (london), 

393, 1443 (1998). 
[11] J.A. Jones and M. Mosca, J. Chem. Phys. 109, 1648 (1998). 
[12] I.L. Chuang, N. Gershenfeld, M. Kubinec, Phys. Rev. Lett. 80, 3408 (1998). 
[13] T. S. Mahesh, Kavita Dorai, Arvind, Anil Kumar, J. Mag. Res. 148, 95 (2001). 
[14] B. M. Fung, Use of pairs of pseudopure states for NMR quantum computing, Phys. Rev. A 

63, 022304 (2001). 

[15] A. K. Khitrin and B. M. Fung, Nuclear magnetic resonance quantum logic gates using 
quadrupolar nuclei, J. Chem. Phys. 112, 6963 (2000). 

[16] T. S. Mahesh, Neeraj Sinha, K. V. Ramanathan, and Anil Kumar, Ensemble quantum- 
information processing by NMR: Implementation of gates and the creation of pseudopure states 
using dipolar coupled spins as qubits , Phys. Rev. A 65, 022312 (2002). 

[17] Neeraj Sinha, T. S. Mahesh, K.V. Ramanathan, and Anil Kumar, Toward quantum informa- 
tion processing by nuclear magnetic resonance: Pseudopure states and logical operations using 
selective pulses on an oriented spin 3/2 nucleus, J. Chem. Phys. 114, 4415 (2001). 

[18] A. Khitrin, H. Sun, and B. M. Fung, Method of multifrequency excitation for creating pseu- 



16 



dopure states for NMR quantum computing, Phys. Rev. A 63, 020301 (2001). 
[19] V. L. Ermakov and B. M. Fung, Experimental realization of a continuous version of the Grover 

algorithm, Phys. Rev. A 66, 042310 (2002). 
[20] A. K. Khitrin and B. M. Fung, NMR simulation of an eight-state quantum system, Phys. Rev. 

A 64, 032306 (2001). 

[21] K.V.R.M. Murali, Neeraj Sinha, T.S. Mahesh, Malcom Levitt, K.V. Ramanathan, and Anil 
Kumar, Quantum-information processing by nuclear magnetic resonance: Experimental im- 
plementation of half- adder and subtractor operations using an oriented spin-7/2 system, Phys. 
Rev. A 66, 022313 (2002). 

[22] Anil Kumar, K.V. Ramanathan, T.S. Mahesh, Neeraj Sinha and K.V.R.M. Murali, Develop- 
ments in quantum information processing by NMT: Use of quadrupolar and dipolar couplings, 
Pramana- Journal of Physics, 59 (2) 243-254 (2002). 

[23] Ranabir Das and Anil Kumar, Use of quadrupolar nuclei for quantum-information processing 
by nuclear magnetic resonance: Implementation of a quantum algorithm , Phys. Rev. A 68, 
032304 (2003). 

[24] R.S. Sarthour, E.R. deAzevedo, F.A. Bonk, E. L. G. Vidoto, T. J. Bonagamba, A. P. 

Guimaraes, J.C.C. Freitas, I.S. Oliveira, Relaxation of coherent states in a two-qubit NMR 

quadrupole system, Phys. Rev. A 68, 022311 (2003). 
[25] F.A. Bonk, R.S. Sarthour, E.R. deAzevedo, J.D. Bulnes, G. L. Mantovani, J.C.C. Freitas, T.J. 

Bonagamba, A. P. Guimaraes, I.S. Oliveira, Quantum-state tomography for quadrupole nuclei 

and its application on a two-qubit system Phys. Rev. A 69, 042322 (2004). 
[26] F.A. Bonk, E.R. deAzevedo, R.S. Sarthour, J.D. Bulnes, J.C.C. Freitas, A.P. Guimaraes, 

I.S. Oliveira, T.J. Bonagamba, Quantum logical operations for spin 3/2 quadrupolar nuclei 

monitored by quantum state tomography, J. Magn. Res. 175, 226 (2005). 
[27] Ranabir Das and Anil Kumar, Experimental implementation of a quantum algorithm in a 

multiqubit NMR system formed by an oriented spin-7/2 system, Appl. Phy. Lett. 89, 024107 



17 



(2006). 

[28] H. Kampermann and W. S. Veeman, Characterization of quantum algorithms by quantum pro- 
cess tomography using quadrupolar spins in solid-state nuclear magnetic resonance, J. Chem. 
Phys. 122, 214108 (2005). 

[29] R. R. Ernst, G. Bodenhausen, and A. Wokaun, Principles of Nuclear Magnetic Resonance in 
One and Two Dimensions, Oxford University Press 1987. 

[30] P. Dhiel and C. L. Khetrapal NMR-Basic principles and progress (Springer- Verlag), New york, 
1969, Vol. 1. 

[31] Y. Aharnov and J. Anandhan, Phase change during a cyclic quantum evolution, Phys. Rev. 
Lett. 58, 1593 (1987). 

[32] D. Suter, K. T. Mueller and A. Pines, Aharnov- Anandhan quantum phase by NMR interfer- 

ometry, Phys.Rev.Lett. 60, 1218 (1988). 
[33] Ranabir Das, S. K. Karthick Kumar, and Anil Kumar, Use of non-adiabatic geometric phase 

for quantum computing by nuclear magnetic resonance, J. Magn. Reson. 177, 318 (2005). 
[34] T. Gopinath and Anil Kumar, Geometric quantum computation using fictitious spin- (1/2) 

subspaces of strongly dipolar coupled nuclear spins, Phys. Rev. A 73, 022326 (2006). 
[35] R. Cleve, A. Ekert, C. Macchiavello, and M. Mosca, Quantum algorithms revisited, Proc. R. 

Soc. Lond. A 454, 339 (1998). 
[36] David Collins, K. W. Kim, and W. C. Holton, Deutsch-Jozsa algorithm as a test of quantum 

computation, Phys. Rev. A 58, 1633 (1998). 
[37] Oliver Mangold, Andreas Heidebrecht, and Michael Mehring, NMR tomography of the three- 

qubit Deutsch-Jozsa algorithm, Phys. Rev. A 70, 042307 (2004). 



18 



(A) 



(B) 



1 7/2, -7/2> 
1 7/2, -5/2> 



t 



g 



1 7/2, -3/2 > 



1 7/2, -1/2 > "■ 

1 7/2, 1/2 > 
1 7/2, 3/2 > 
1 7/2, 5/2 > 
17/2, 7/2 > 



■ ■ 

't 

■ ■■■ ' 

!5 



u 



u 



UL 



— i — 

2e+04 



1 

Oe+00 



Hz 



I" 



|1H>= |8> 
II 1 0> = |7> 
|101>= |6> 
|100>= |5> 
|011> = |4> 
|010>= |3> 
10 1> = |2> 
|000>= |1> 

FIG. 1: (A) Schematic energy level diagram of a spin-7/2 nucleus, spin states 1 7/2, 7/2) 

1 7/2, —7/2) are labeled as eight basis states of a three qubit system, the dark rectangles represent 

the equilibrium populations, and the single quantum transitions are labeled as a, b...., g. (B) Equi- 
librium spectrum of oriented 133 Cs nucleus (spin-7/2), obtained by a non-selective (vr/2) y pulse, 
at a resonance frequency of 65.59 MHz on a Bruker AV-500 spectrometer at a temperature 307K. 
The distance between successive transitions is equal to the effective quadrupolar coupling (A). The 
line widths of transitions a, b, ...,g are observed to be in the ratio 3.2:2.3:1.5:1.0:1.5:2.2:3.1 and 
the integrated experimental intensities are in the ratio 7.2:12.1:15.1:16:15.1:12.0:7.0 (theoretically 
expected 7:12:15:16:15:12:7) 
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FIG. 2: A two level subspace (r,s) is represented on a Bloch sphere. Two it pulses applied on a 
transition (r,s) with phases 6 and (9 + tt + </>), cyclically rotates the states |r) and |s). Solid angle 
subtended by this closed loop at the center of the sphere is, 2<p. The phase acquired by the states 
\r) and \s) is e 1 ^ and e~ l< ^ respectively. 
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FIG. 3: Implementation of controlled phase shift gates (a) UL ^ = diag[l, 1,-1,-1,1,1,1,1] and 
(b) U^ 2 = diag[e t7r ^ 2 , 1, 1, 1, 1, 1, 1, 1], preceded by a hard (vr/2)^ pulse. is implemented by 

applying two (ir) y pulses on transition c. U^ 2 is implemented by applying four multi-frequency it 
pulses (Eq. BJ. 
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FIG. 4: Schematic energy level diagram of an N-level system. The arrows represent single quantum 
transitions. 
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FIG. 5: Quantum circuit for Collins version of DJ algorithm on a three qubit system. 
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FIG. 6: (a) Population distribution of 1 000) PPS obtained by applying an amplitude modulated 
multi-frequency pulse (explained in text) on transitions b, c, d, e, f, and g at equilibrium state, 
(b) Spectrum of |000) PPS obtained by applying a non selective 5° pulse with 32 scans, (c) 
Coherent superposition state (Eq.s [7|8]) obtained by applying a (ir/2) y pulse on 1 000) PPS. The 
intensities are in accordance with Eq. ([7]), which is different from equilibrium intensities of Fig. 
[TJ Experimental integrated intensities are in the ratio 10:48:101:140:100:46:8, while the expected 
theoretical intensities are in the ratio 7:42:105:140:105:42:7. The intensities here represent unequal 
superposition of basis states. 
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FIG. 7: Implementation of Collins version of 3-qubit DJ algorithm (Fig. [5]). Spectrum oiipd which 
corresponds to a constant function U=I. Spectra of ipikjm (Eq. [12]) which correspond to balanced 
functions U(l,k,l,m) (Eq. [6]). The required MF pulses for implementing U(l,k,l,m) are given in 
Table (2). The duration of each MF pulse is 1.425ms. Each spectrum is recorded in 4 scans 



